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The intersection volume of two independent 2-level cut 
Gaussian random fields is proposed to model the open-cell 
microstructure of organic aerogels. The experimentally mea- 
sured X-ray scattering intensity, surface area and solid ther- 
mal conductivity of both polymeric and colloidal organic aero- 
gels can be accounted for by the model. 
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Aerogels are a promising material for a host of applica- 
tions [|l],^| due to their thermal, optical and mechanical 
properties. For example, aerogels are among the best 
thermal insulating solid materials known (^-(^]. It is im- 
portant to link aerogel properties to their complex inter- 
nal microstructure, and to understand how such proper- 
ties can be optimized for a given application [HUH • The 
nano-scale porous morphology of aerogels has been ex- 
tensively characterized by X-ray scattering and surface 
area analysis [|7|— [Toj] . Despite this, aerogel properties are 
usually correlated with density, rather than related to 
morphological features. One reason for this is the lack 
of suitable representation of aerogel morphology. In this 
paper we develop a model which accounts for the open- 
cell morphology of organic aerogels. The solid thermal 
conductivity of the model is computed and shown to be 
in good agreement with experimental data. 

Thermal transport in aerogels is due to three additive 
components: conduction in the solid skeleton and (gas- 
filled) pores and conduction due to radiation j^]. For 
thermal insulation purposes it is desirable to reduce the 
magnitude of each contribution. Gaseous conductivity 
can be significantly reduced by decreasing pore size or 
partially evacuating the material, and radiative transport 
reduced by the inclusion of an opacifier (^^]. The solid 
conductivity (typically half the total) depends strongly 
on the aerogel density and microstructure . 

Organic aerogels produced by the polymerisation of re- 
sorcinal and formaldehyde (RF) have been suggested as 
an alternative insulator to opacified silica aerogels ||,[ll| . 
They have lower intrinsic and radiative conductivities, 
and are less brittle than their silica based counter- 
parts J|,|]]. Both the morphology [j3||T^1 and properties 
of organic aerogels have been the subject of 

detailed investigation. A key variable in the formation 
of RF aerogel microstructure is the initial ratio of rc- 
sorcinal to catalyst (R/C) ||. As the catalyst increases 
the aerogels vary from a colloidal structure to a well- 
connected polymeric structure with a corresponding in- 



crease in conductivity and strength [|j,|13| . It is important 
to quantitatively model these properties to assist in the 
understanding and optimization of RF aerogels. 

Current models of aerogels are based on simulating 
the microstructure formation using the diffusion-limited 
cluster-cluster aggregation (DLCA) scheme p"8| . 
Two features of DLCA models (proposed for silica aero- 
gels) suggest that they are not well suited to modeling 
RF aerogels. Firstly, the DLCA model exhibits fractal 
scaling Jl(|[l7j and a well pronounced peak in the scat- 
tering intensity jlO 18|. In contrast, RF aerogels exhibit 
no fractal scaling, and under high catalyst conditions the 
peak is weak, or even absent ||. Secondly, the discrete 
character of DCLA type models (open networks of cubes 
or hard spheres) may be ill-suited to modeling "contin- 
uum" properties within the aerogel skeleton. For exam- 
ple, the influential inter-particle neck size fllij is equal 
to zero for hard spheres and equal to the particle 
size for cubes ]l6|-p^|. We propose a statistical model 
of microstructure which can account for the main mor- 
phological features of RF aerogels. The model is lattice 
independent, and suitable for continuum based theoreti- 
cal and computational prediction of properties jl9|,^0| . 

A convenient statistical description of porous media is 
provided by modeling the internal interface as the iso- 
surface (or level-cut) of a Gaussian random field (GRF) 
y(r). This model has been used to describe the mor- 
phologies arising in spinodal decomposition |2~l[ , mi- 
croemulsions p2| , p3| , galaxy formation p4| and porous 
rocks |£5| amongst others [^6| . The statistics of the mate- 
rial are completely determined by the specification of the 
single level-cut parameter and the field-field correlation 
function g(r) — (y(O)y(r)) (where r — \r\ and <?(0) = 1). 
Berk generalized the model to account for the X-ray scat- 
tering properties of microemulsions p7| by defining phase 
1 to occupy the region in space where a < y(r) < (3 and 
phase 2 to occupy the remainder. This model has also 
been shown to account for the morphology and properties 
of foamed solids |^8| and polymer blends . 

Neither the 1-cut GRF model, or Berk's "2-cut" ex- 
tension, can account for the high porosity open-cell mi- 
crostructure of aerogels. The 1-cut GRF is not macro- 
scopically connected at aerogel porosities Q (typi- 
cally 95%), and Berk's 2-cut model exhibits sheet-like 
structures |3(| similar to those observed in closed-cell 
foams To model the open-cell microstructure we 

define the solid phase to occupy the region a < y(r) < (3 
and a < w(r) < (3 where y and w are statistically inde- 
pendent GRF's. The independence of the random fields 
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allows the correlation functions of the model to be cal- 
culated. The solid volume fraction of the model is just 
p = (pp — pa) 2 where p a — (2-k)~^ J" e - * l 2 dt and the 
2-point correlation function is P2 (r) = h\ (r) where hi (r) 
is the 2-point function of Berk derived for the same values 
of a and (3 
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with = pp —p a . The freedom in specifying the level-cut 
parameters and the field-field correlation function of the 
model g(r) allow a wide variety of morphologies to be 
modelled. 

To relate the model to experimental data it is neces- 
sary to specify a field-field correlation function. Prior 
studies p2|32j suggest a form 
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characterized by a correlation length £, domain scale <i 
and a cut-off scale r c . Two commonly measured morpho- 
logical quantities of porous media are the surface area 5* 
and the X-ray scattering intensity I(q). These can be 
computed for the model as 
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where V is the sample volume and rj is the scattering 
density of the solid phase ||| . 

To model RF aerogels we choose the model parameters 
to match experimentally measured scattering and surface 
area data B. While the domain scale d corresponds to 
the pore scale in aerogels, the geometry of the fibres de- 
pends on both the length scale and level-cut parameters. 
Uncertainties in the estimation of surface area ||[l2],[34| 
and skeletal density of aerogels suggest that only rough 
approximations of the parameters are justified. Examples 
of a colloidal and polymeric aerogel are chosen to ascer- 
tain the generality of the model. The colloidal aerogel is 
produced under low catalyst conditions (R/C=300) and 
has density p a = 148kg m - 3 (p = p a /p s = 0.11) §. TEM 
images show a "string-of-pearls" appearance: a network 
comprised of grains (diameter 120-300A) interconnected 
by narrow necks. The surface area is 400m 2 /g and X-ray 
scattering yields a peak at q = 0.012A -1 (Fig. @) @. A 
good match between the experimental information and 
the model was obtained for r c — 10, £ = 14 & d — 46nm 



and p a .fi = .07, .40. The model surface area is 428m 2 /g, 
and the theoretical scattering curve is seen to be in good 
agreement with the experimental data (Fig. |l|) . In Fig. || 
we show a slab of the material: the model reproduces 
the colloidal string-of-pearls morphology with reasonable 
grain and neck sizes. We directly measure the scatter- 
ing from the model. The results are included in Fig. [l], 
demonstrating that the simulation reproduces the theo- 
retical statistics very well. At small q (large length scales) 
some deviation is evident; this is due to the finite size of 
the samples p5| . 

The polymeric aerogel || we model is produced un- 
der high catalyst conditions (R/C=50) and has density 
p a = 100kg/m 3 (p — 0.077). The aerogel exhibits a net- 
work of uniform fibres (diameter 30-60A) with surface 
area 905m 2 /g. The scattering intensity monotonically 
decreases with q and yields a Guinier radius of 81 A [||. 
In this case we take r c = 10,£ = 20&<i = 30nm and cen- 
tered level cut parameters (a — —f3) p a ,(3 — .361,-639. 
The model has a surface area of 927m 2 /g and a Guinier 
radius of lOlA. Fig. [l] shows good agreement between 
the model and experimental scattering curves. A slab of 
the model material is shown in Fig. ^. The fibres have a 
relatively uniform thickness, varying from 20 — 60A. 

As the presence of a peak in the scattering may yield 
information about the the physical processes underlying 
aerogel formation p|fL^] it is interesting to comment on 
its morphological origins. The existence of a domain (or 
repeat) scale in a random structure leads to decaying 
oscillations in the correlation function, and hence a peak 
in I(q). In aerogels the decay scale is controlled by the 
width of the fibres Wf, and the domain scale d is that 
of the pores. If d is only several times larger than wj 
(e.g. the colloidal model) a peak is observed. On the 
other hand, if d is an order of magnitude larger than 
Wf (as it is in the polymeric model) the oscillations in 
Pi{r) are smoothed by a stronger decay and the peak is 
extinguished (Fig. [l]). Note that pores with a well defined 
scale are evident in the model (Fig. |J); they simply do 
not carry sufficient statistical weight to appear in the 
scattering. 

We now compare the thermal conductivity of the 
model to experimental data. The solid thermal conduc- 
tivity of RF gels has been experimentally measured over 
the density range p a = 60 — 400kgm -3 at catalyst con- 
centrations (R/C=50,200 & 300) |[§]. To estimate the 
model conductivity we assume that the local heat flux 
obeys the Fourier law j — — A(r)VT where A(r) = A s (0) 
in the solid(void) phase. Conventional numerical tech- 
niques are used to solve the heat conservation equation 
V • (AVT) = in a 128 3 lattice subject to an applied 
temperature gradient (^0|. The aerogel conductivity is 
obtained as A a = (A(r)VT)/(VT). 

Note that we have derived models of colloidal and poly- 
meric aerogels based on the morphology and scattering 
data at a specific density. We can extend the model to 
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higher and lower densities by making simple assumptions 
about the density dependence of the model parameters. 
A simple scaling argument shows the pore scale varies 
as d oc p~z pj and for simplicity a similar dependence 
is assumed for r c and £ (the conductivity is relatively 
insensitive to length scale variations |^,|o|). The poly- 
meric morphology of the high catalyst (R/C=50) aerogels 
was reproduced by "centering" the level cut parameters 
(a = —j3) (Fig. ||). We preserve this feature of the model 
by choosing pp — \ + hs/P and p a = \ — \ \[V- The ther- 
mal conductivity of the polymeric model "P" is shown 
in Fig. ^. While slightly underestimating the conductiv- 
ity of the aerogel produced at R/C=50, it is nevertheless 
seen to be in very good agreement with the experimental 
data. At low catalyst concentration (R/C=300) the aero- 
gel morphology was modeled by non-centered level cuts 
(a 7^ —jS) (Fig. ^|). To maintain an asymmetry we choose 
Pfj = 0.3+ ^y/p and p a = 0.3— ^^/p. The thermal conduc- 
tivity of this model "C" is presented in Fig. ^, providing 
excellent agreement with the experimental data for the 
RF aerogel produced at R/C=300. 

In Fig. U we have also plotted a number of results 
arising from theoretical considerations. Zeng et al |36| ] 
have suggested that periodic open-cell models can be 
used to estimate aerogel conductivity. At low relative 
densities the "square rod" model leads to the estimate 
A Q /A S = \pal ' Ps ||o) in remarkably good agreement with 
the data for polymeric RF aerogels. From considerations 
of phonon heat transport in solids it has been suggested 
that A a /A s = p a Val PsV s where v a (v s ) is the sound veloc- 
ity in the aerogel(solid) J37J. Measurements performed on 
RF aerogels have determined that v a /v s — 0.47(p/p s ) ' 88 
so that Aa/% = 0A7(p a /p s ) im f§[||. The result is 
seen to significantly under-estimate the measured con- 
ductivity of RF aerogels. It is also possible to calculate 
rigorous variational bounds j3j],^(J on the model ther- 
mal conductivity |20|j3C| : |35| . The upper bounds, which 
can have predictive power [|l9U3C| 1, are seen to consider- 
ably over-estimate the true conductivity. Thus theoret- 
ical microstructure-property relations are unable to pre- 
dict the thermal conductivity, and numerical simulations 
must be relied on. 

The agreement between the model and experimental 
data for both colloidal and polymeric aerogels provides 
strong evidence that we have accurately modelled the 
morphology of organic aerogels. The results also indicate 
that Fourier's continuum theory of heat conduction may 
hold even in nano-scale structures (diameter 30-60 A). Of 
course there is no guarantee that the model is correct: 
other models may share the same morphological pjq] and 
thermal properties. Nevertheless the utility of the model 
has been shown. It should be possible to apply the model 
in the study of gas and radiative conductivity and the 
mechanical properties of aerogels. The fractal properties 
of silica aerogels can also be incorporated |35| . 

Extensions of the model are relevant to a wider range 



of heterogeneous materials. For example, the solid phase 
of the aerogel model mimics the inter-granular pores of 
sandstone, and micro-porosity can be simulated by in- 
cluding random structures at smaller scales. Spheres 
may be embedded in the models, and closed cell mor- 
phologies, such as those observed in solid foams can 
be formed from the union set of two 2-level cut GRFs. 
The model correlation functions can be calculated, allow- 
ing surface areas, scattering curves and rigorous property 
bounds to be evaluated [05|. 
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FIG. 1. X-ray scattering spectra of RF aerogels. The data 
for the upper and lower curves are for a colloidal (R/C=300, 
p=ll%) and polymeric (R/C=50, p—7.7%) aerogel respec- 
tively |8| (the experimental data are vertically scaled). The 
theoretical curves were obtained using Eq. (Q) , and the simu- 
lation data are measured directly from one realization of each 
model (Figs. | & |). The inset shows log 10 {I(q)/V(T] 2 )) vs. 
!ogi (9)- 




FIG. 2. A model colloidal aerogel (solid fraction 
p =11%). The parameters are r c =10, £=14, d=46nm and 
p a ,f3 = 0.07,0.40. The image is 276 x 276 x 34.5nm. The 
slab is part of a periodic cubic sample of side length 276nm 
(128 pixels). Many of the apparently isolated clusters are 
interconnected outside the volume shown. 




FIG. 3. A model polymeric aerogel (solid fraction 
p — 7.7%). The model parameters are r c =10, £=20, d=30nm 
and p a ,0 = 0.36, 0.64. The image is 90 x 90 x 11.25nm. 
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FIG. 4. Solid thermal conductivity of RF aerogels; theory 
vs. experiment g] (p s = 1300kg/m 3 , X a = O.lSWm^KT 1 }. 
The estimates of A a from this work (solid symbols) show 
very good agreement with experimental data (open symbols). 
The solid and dashed lines correspond to predictions of X a 
discussed in the text. The dotted lines are rigorous upper 
bounds. 
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